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ABSTRACT 

We consider the problem of estimating the virial mass of a dark halo from the po- 
sitions and velocities of a tracer population. Although a number of general tools are 
available, more progress can be made if we are able to specify the functional form of 
the halo potential (although not its normalization). Here, we consider the Navarro- 
Frenk- White (NFW) halo and develop two simple estimators. We demonstrate their 
effectiveness against numerical simulations and use them to provide new mass esti- 
mates of Carina, Fornax, Sculptor, and Sextans dSphs. 

Subject headings: galaxies: halos — galaxies: kinematics and dynamics — dark mat- 
ter 

1. Introduction 

There are many instances when we wish to estimate the mass of a dark halo based on proper- 
ties of a tracer population. This includes calculating the masses, of the Milky Way and M3 1 halos 
from the positions and velocities of satellite galaxies, globular clusters and halo stars; of the dark 
halos of dwarf spheroidal galaxies from the projected positions and velocities of constituent stars; 
of galaxy clusters from the motions of member galaxies. 

There has be en much previous work on this problem. Early estimators were based on the 



virial theorem, but Bahcall & Tremaine (1981) showed that such estimators are inefficient, biased 



and, inconsistent, and replaced them with the projected mass estimator. Their original estimator 
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was tailored for test particles around a poi nt mass whereas IHeisler et al.l (|1985|) extended it for 
self-consistent cases. This was modified by IWatkins et al.l (|2010f) to deal with tracer populations 
whose number dens ity does not necessarily trace an extended underlying dark matter distribution. 
An & Evansl (1201 1|) on the other hand devised a general theory of mass estimators under the as- 
sumption that the dark matter density law is known, although the normalization is not. 

Here, we start by a ccepting that numero us studies have shown us that the dark halos have the 
universal form, namely (|Navarro et al.lll995i henceforth NFW), 
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where a is the scale radius. We wish to find the virial mass M yn , which is the mass within the virial 
radius r vir . We deduce the enclosed mass M{r) within the radius r; 
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and the normalized mass profile M(r) follows as 
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Here M out = M(r out ), and r out is the outermost data point in our sample while f = r/a = cr/r vir , 
and c = r vir /a is the concentration parameter. We have here assumed r out = r vir and so M out = M v i r , 
which is the case for halo mass estimation using satellite galaxies. 

The problem therefore is: Given an assumed NFW profile together with the positions and 
velocities of tracers that extend out to the virial radius, how can we estimate the virial mass? 
The aim of this Letter is to provide a ready-to-use specifically-optimized estimator for the most 
popular halo model with tests and applications, without burdening the reader with a plethora of 
mathematical derivations. 



2. Two Estimators 



An & Evansl (120111) show that if we are willing to assume the functional form of the halo mass 
profile M(r) within the given spherical region of interest, then the total dark halo mass M vir in the 
same region can be estimated through a weighted average of kinematic properties of the tracers, 
adjusted by the boundary term. The boundary term is related to the external pressure support of the 
tracer population, but it is usually taken as small compared to other uncertainties and so neglected. 
The mass of the NFW halo is formally infinite. This means that the profile needs to be truncated 
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at some radius, which is usually taken to be either r out or r vlT . In this paper, we suppose that the 
tracers are well-populated out to the virial radius so that r out = r vir . 

We consider two cases. First, as for the Milky Way satellites, 3-d distances r to and radial 
velocities v r with respect to the halo center are available for tracers. Strictly speaking, the radial 
velocity is with respect to the Sun, but for most practical proposes, this is almost equivalent to the 
the radial velocity with respect to the Galactic Center for distant tracers. Second, as for populations 
in most external galaxies, projected distances R and line-of-sight velocities v c are available for 
tracers. We consider these two in turn. 



2.1. True Distances and Radial Velocities 



The general principle underlying all estimators is that the mass is a weighted average of the 
positions and velocitie s of th e tracers. For a given choice of functional form for the dark matter 
potential, lAn & Evansl (|201l|) showed that there is an optimum weight function. For the case of 
this NFW halo, this gives 
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Here, J3 is the Binney anisotropy parameter for the spherical system. Deason et alJ (|2011r) analyzed 
properties of velocity distributions of satellites galaxies in numerical simulations of dark halos and 
found that typically ft « 0. We retain f3 in our formulae, but in applications we assume (3 = 0. 
The last term in equation © is the boundary condition, which is discarded here. This gives the 
estimator as 
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where the index i runs over the Af tracers. We refer to this as the NFW Mass Estimator (NFWME). 



It is useful to compare this with a special case of the estimator derived by IWatkins et al. 



They studied dark halos with density profiles that are scale-free out to the virial radius 



/ r \ l ~ a 
M(r)= - • 

\ / vir ' 

For this case, the weight function is a power-law, and the mass estimator becomes 

(a + 7 -2j3)r l - a 
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where y is the power-law index for the tracer density (oc r y ). [An & Evan s (201 1) have shown tha t 



the choice y = 3 corresponds to neglecting the boundary term. By contrast, IWatkins et al.l (120 lQf) 
circumvented the effect of the boundary term by explicitly solving for the power- law solution for 
the tracers. Here, we recommend the choice y « 2.5 for satellite galaxy populations, but y « 3 for 
stellar populations in dark haloes. Simulations ty pically find that th e number densit y of satellite 



galax ies falls off more slowly than the stars (e.g.. iDeason et al.l 1201 II) . In addition. IWatkins et al. 
(1201 Oh also argued that a ~ 0.5 is a good approximation to the NFW halo, at least for the range 
of virial masses and concentrations appropriate to large galaxies like the Milky Way. With these 
particular indices, the mass estimator is in the form: 
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We refer to this as the Scale-Free Mass Estimator (SFME). Equations © and © are the funda- 
mental results of this subsection. They provide simple formulae for the mass in terms of weighted 
sums of the positions and velocities. We shall test these formulae against simulations shortly. 



2.2. Projected Distances and Line-of-Sight Velocities 

In many situations, the radial velocities and true positions of the tracers with respect to the 
center of the halo are not direct observables, but the line-of-sight velocities and projected posi 



tions are. The adjustment of the estimator for this happenstance is covered in lAn & Evansl (|201 1 
Sect. 4). 

Here, we would like to find the weight function W(R) of R such that (y iW(R)) = {rw(r)M~ l v 2 r ), 



which is to be substituted into the mass estimator. lAn & Evansl (|2011l) show that this leads to 
an integral equation and the weight function W(R) can be computed, at least numerically. The 
invariant (i.e., independent of c) normalized weight functions W(R) = W(R)/[a^(c)] for the NFW 
profile for some constant ft are shown in Figure Q] 

Notice that for much of the range occupied by the tracers and of the physical range of the 
anisotropy parameter (l/2>/3>-lor the axis ratio of the velocity ellipsoid no more extreme 
than 1 : V2), the weight function in Figure [T]is constant to a good approximation, W(R) « 15, so 
that 

M ( C ) GN 4-i K a> ° GN 4^ ° 

(=i i=i 

Interestingly, this indicates that M(a) « 3G~ l a{v 2 i ) since M vir //i(c) = M(r)//u(r/a) and yu(l) « 0.2. 
In what follows, we however use the numerically obtained weight function for the mass estimator, 
which is the projected analog of equation ©. 
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We will again compare this estimator to the projected analog of the scale-free estimator. If /3 
is a finite constant a nd the mass profile is in the scale-free form in equation ©, we find that (c.f., 



Watkins et alJbOld eqs. 26 & 27) 
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Specializing to the case a « 0.5 and y « 3, we find 
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This is the projected analog of equation ©. Equations (flOl and (fT2j) are the fundamental results 
of this subsection. 



3. Applications 



3.1. Numerical Simulations 



We begin by testing our estimators against simulations. The Galaxies-In tergalactic Medium 
Interaction Calculation (GIMIC) suite of simulations is described in detail in ICrain et al.l (|2009|) . 
It consists of a set of hydrodynamical resimulations of five nearly spherical r egions (~ 20/T 1 Mpc 
in radius) extracted from the Millennium Simulation ([Springel et al. 2005 ). Deason et al. ( 2011 ) 
extracted from the GIMIC simulations a set of galaxies that resemble the Milky Way. The sample 
consists of 431 parent halos and 4,864 associated satellite galaxies. 

There are a number of ways in which mock catalogs from simulations differ from the as- 
sumptions used to derive the estimators. For example, dark halos are not generally spherical, 
infall continues to the present day, and the observed satellites are not necessarily virialized and 
well-described by an equilibrium distribution. Also, notions of a constant anisotropy are probably 
idealized, and in practice the anisotropy will vary with radius. These systematic errors are larger 
than random errors on the measurements of velocities and positions. 

The mass estimators provide an estimate for the total mass within the radius of the farthest 
tracer (« r out ). We compute the 'true' mass within r out « r vir for each halo and compare to masses 
found via our estimators. We use all satellites, but check that our results are not significantly 
affected when only luminous satellites are included. Results are obtained for the two estimators 
in this paper - the NFWME and th e SFME - and also t he vir ial mass estimator (VME) and the 



projected mass estimator (PME) of iBahcall & Tremainel (|198l|) . The left panel of Figure prefers 
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to data sets of radial velocities and true distances, and the right panel to line-of-sight velocities 
and projected distances. We quantify the goodness of any estimator by means of two statistical 
measures. First, we define the Fraction of Reasonabl e Estimates (FRE) a s the fraction of estimates 
within the factor of two of the true mass (see also Deason et al.ll201 We also give the Inter 
Quartile Range (IQR) of the mass estimates, which gives a good indication of the spread. In 
all cases, we see that the NFWME and the SFME outperform the VME and PME. Interestingly, 
the performance of the NFWME is comparable to SFME. Given this, we recommend the use of 
equations © and (PT21) for practical applications. 



3.2. Dwarf Spheroidal Masses 



We now turn to an astrophysical problem. The recent years have seen programs to garner the 
radial velocities of giant stars in the nearby dwarf spheroidal galaxies (dSphs). For the brighter of 
the dSphs, radial velocity surveys have pr ovided data sets of projected positions and line-of-sight 



velocit ies for thousands of stars (see e.g.. Kleyna et al.ll2002 ; Wilkinso n et al.l 12 004: 



Walker et al. 



2009a ). This has driven new theoretical ideas, such as the claim by lstrigari etaL (|2008 ) that all the 
dSphs shared a common mass scale of ~ 10 7 M Q with in 300 pc. The data has also driven the stud y 



of new techniques for mass estimation and modeling (|Wolf et al.ll2010l : lAmorisco & Evansll201 1|) . 



We use the velocities and positions of individu al stars observed in four classical dSphs - 
Carina, Fornax, Sculptor, and Sextans - presented by IWalker et al.1 (I2009a|) . We select only those 
stars for which the probability of membership is > 0.9, and additionally assume that the isotropy 
(fi = 0). The results for the mass within 300 and 600 pc are listed in Table [D There are a couple of 
points worth noting. First, the masses do indeed support lStrigari et al.l(|2008|) 's notion of a common 
mass scale within 300 pc, although this is less surprising as the velocity dispersions of these four 
dSphs are similar. Second, the PME agrees well with the two new estimators for the enclosed mass 
within 300 pc. However, it begins to diverge from the SFME and NFWME at larger radii. 

Although the re sults in Table \T\ are compa rable to mass estimates obtained through Jeans 
(IWalker et al.ll20 09b) and distribution function (|Amorisco & Evansll2011r) modeling, the amount 
of effort involved is very much less. All that is required is a weighted average of positions and 
velocities as opposed to solving differential, or integro-differential equations. 



4. Conclusions 



We have provided a new and accurate way of estimating the masses of NFW halos from the 
positions and radial velocities of tracer populations. This work follows up the theoretical paper of 



An & Evansl (|2011|) and supplies a ready-to-go mass estimator for the most common application 
of all. 



We have supposed that the reader has a data set of N objects with true positions r, and radial 
velocities v r j, or with projected positions i? ( and line-of-sight velocities v C j. In terms of overall 
simplicity and flexibility, we recommend using the isotropic limit (J3 = 0) with an r" 3 tracer number 
density fall-off (y = 3) 

]U ~ J - jr out Y 0.5 2 

' =1 (13) 

4 0.5 n 

!=1 

to compute the enclosed mass within the sphere of the radius r out . Here, r out is the location of the 
outermost data point. In the projected case, only R om is known, but on statistical grounds, we have 
that r out = nR out /4. Obviously, no information can be inferred on the mass distribution exterior to 
the outermost data point. 
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Fig. 1. — The invariant normalized weight function W(R) for the NFW profile as a function of the 
normalized projected separation R = R/a. The weight function can be series-expanded at R = 
with its leading term behaving as ~ R~ l unless f3 = 1. On the other hand, asymptotically toward 
R -» oo, we find that W(R) ~ R/ log J?. 
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Fig. 2. — Histograms showing the ratio of the estimated mass to the true mass for the two different 
mass estimators. The left panel uses eqs. © and ® applied to data sets of true distances and radial 
velocities. The right panel uses eqs. (flOl) and (fT2l) applied to data sets of projected distances and 
los velocities. For the NFWME estimator, a concentration c = 17 is assumed, which is typical for 
the haloes in the simulations. FRE and IQR are the fraction of reasonable estimates (FRE) and 
the inter-quartile range (IQR) as defined in the text. The means of the distributions are shown by 
vertical dotted lines, from which we see th at the virial mass estimator is biased and underestimates 
the true mass (as originally pointed out by iBahcall & Tremaindll98l|) . 
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Table 1. Dwarf Spheroidal Masses 





dSph 


PME 


SFME NFWME 

(X 10 7 Me) 


Carina 


(R < 300 pc) 


1.7 


1.7 


1.7 




(R < 600 pc) 


2.2 


2.6 


2.7 


Fornax 


(R < 300 pc) 


2.9 


2.8 


2.9 




(R < 600 pc) 


5.1 


5.2 


5.2 


Sculptor 


(R < 300 pc) 


1.8 


1.8 


1.8 




(R < 600 pc) 


2.5 


3.1 


3.0 


Sextans 


(R < 300 pc) 


1.3 


1.3 


1.3 




(R < 600 pc) 


2.0 


2.2 


2.2 



Note. — All the error are +0.1 x 10 7 M Q , taking into 
account the random measurement uncertainties only. 
For the NFWME, a concentration c = 30 (typical for 
dSphs) has been assumed. 



